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Abstract 

A unifying approach to software and hardware design generated by 
ideas of Idempotent Mathematics is discussed. The so-called idempo- 
tent correspondence principle for algorithms, programs and hardware 
units is described. A software project based on this approach is pre- 
sented. 
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1 Introduction 

Numerical computations are still very important in computer applications. 
But until recently there was a discrepancy between numerical methods and 
software/ hardware tools for scientific calculations. In particular, numeri- 
cal programming was not much influenced by the progress in Mathematics, 
programming languages and technology. Modern tools for numerical calcu- 
lations are not unified, standardized and reliable enough. It is difficult to 
ensure the necessary accuracy and safety of calculations without loss of the 
efficiency and speed of data processing. It is difficult to get correct and exact 
estimations of calculation errors. For example, standard methods of interval 
arithmetic \W\ do not allow to take into account the error auto-correction 



effects [jT9[ and, as a result, to estimate calculation errors accurately. 

However, new ideas in Mathematics and Computer Science lead to a very 
promising approach (initially presented in [^-[Q). An essential aspect of 
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this approach is developing a system of algorithms, utilities and programs 
based on a new mathematical calculus which is called Idempotent Analysis 
or Idempotent Calculus, or Idempotent Mathematics etc. For many problems 
in optimization and mathematical modeling this Idempotent Calculus plays 
the same unifying role as Functional Analysis in Mathematical Physics, see, 
e.g., ig, 0, and surveys 

Idempotent Analysis is based on replacing the usual arithmetic operations 
by a new set of basic operations (such as maximum or minimum). There are 
a lot of such new arithmetics which are associated with sufficiently rich alge- 
braic structures called idempotent semirings. It is very important that many 
problems, nonlinear in the usual sense, become linear with respect to an ap- 
propriate new arithmetic, i.e. linear over a suitable semiring (the so-called 
idempotent superposition principle p6[, [0, fl^, which is a natural analog of 



the well-known superposition principle in Quantum Mechanics). This 'linear- 
ity' considerably simplifies explicit constructions of their solutions. Examples 
are the Bellman equation and its generalizations, the Hamilton-Jacobi equa- 
tion etc. The idempotent analysis is a powerful heuristic tool to construct 
new algorithms and apply unexpected analogies and ideas borrowed, e.g., 
from mathematical physics and quantum mechanics. The abstract theory 
is well advanced and includes, in particular, a new integration theory, lin- 
ear algebra and spectral theory, idempotent functional analysis, idempotent 
Fourier transforms and so on. Its applications include various optimization 
problems such as multi-criteria decision making, optimization on graphs, dis- 
crete optimization with a large parameter (asymptotic problems), optimal 
design of computer systems and computer media, optimal organization of 
parallel data processing, dynamic programming, applications to differential 
equations, numerical analysis, discrete event systems, computer science, dis- 



crete mathematics, mathematical logic, etc. (see, e.g. pi], H], [H]-[|T5|) ¥]_ 



[21], [26]-] 32] and references therein). 

It is possible to obtain an implementation of the new approach to sci- 
entific and numeric calculations in the form of a powerful software system 
based on unified algorithms. This approach ensures the arbitrary necessary 
accuracy and safety of numerical calculations and a working time reduction 
for programmers and users because of a software unification. 

Our approach uses the techniques of object oriented and functional pro- 
gramming (see, for example, [^, [0) which is very convenient for the design 
of our (suggested) software system. A computer algebra techniques M is also 
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used. The modern techniques of systohc processors and VLSI reahzations of 
numerical algorithms including parallel algorithms of linear algebra (see, for 
example, |]TB|], |^) is convenient for effective implementations of the pro- 



posed approach to hardware design. 

There is a regular method based on the theory for constructing back-end 
processors and technical devices intended for a realization of basic algorithms 
of idempotent calculus and mathematics of semirings. These hardware facil- 
ities can increase the speed of data processing. 



2 Mathematical objects and their computer 
representations 

Numerical algorithms are combinations of basic operations. Usually these ba- 
sic operations deal with 'numbers'. In fact these 'numbers' are thought of as 
members of some numerical domains (real numbers, integers etc.). But every 
computer calculation deals with finite models (or finite computer represen- 
tations) of these numerical domains. For example, integers can be modeled 
by integers modulo a power of number 2, real numbers can be represented 
by rational numbers or floating-point numbers etc. Discrepancies between 
mathematical objects (e.g. 'ideal' numbers) and their computer models (rep- 
resentations) lead to calculation errors. 

Due to imprecision of sources of input data in real-world problems, the 
data usually come in the form of confidence intervals or other number sets 
rather than exact quantities. Interval Analysis (see, e.g., 0) extends op- 
erations of traditional calculus from numbers to number intervals to make 
possible processing such imprecise data and controlling rounding errors in 
computational mathematics. 

However, there are no universal models which are good in all cases and 
we have to use varieties of computer models. For example, real numbers can 
be represented by the following computer models: 

standard floating-point numbers, 

double precision floating-point numbers, 

arbitrary precision floating-point numbers, 

rational numbers, 

finite precision rational numbers. 
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floating-slash and fixed-slash rational numbers, 
interval numbers, etc. 

To examine an algorithm it is often useful to have a possibility to change 
computer representations of input/output data. For this aim the correspond- 
ing algorithm (and its software implementation) must be universal enough. 

3 Universal algorithms 

It is very important that many algorithms do not depend on particular models 
of a numerical domain and even on this domain itself. Algorithms of linear 
algebra (matrix multiplication. Gauss elimination etc.) are good examples 
of algorithms of this type. 

Of course, one algorithm may be more universal than another algorithm 
of the same type. For example, numerical integration algorithms based on 
the Gauss- Jacobi quadrature formulas actually depend on computer models 
because they use finite precision constants. On the contrary, the rectangular 
formula and the trapezoid rule do not depend on models and in principle can 
be used even in the case of idempotent integration (see below). 

The so-called object oriented software tools and programming languages 
(hke C^^ and Java, see, e.g., [^) are very convenient for computer imple- 
mentation of universal algorithms. 

In fact there are no reasons to restrict ourselves with numerical domains 
only. Actually it may be a ring of polynomials, a field of rational functions, 
or an idempotent semiring. The case of idempotent semirings is extremely 
important because of numerous applications. 

4 Idempotent correspondence principle 

There is a nontrivial analogy between Mathematics of semirings and Quan- 
tum Mechanics. For example, the field of real numbers can be treated as 
a 'quantum object' with respect to idempotent semirings. So idempotent 
semirings can be treated as 'classical' or 'semi-classical' objects with respect 
to the field of real numbers. 

Let R be the field of real numbers and R+ the subset of all non-negative 
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numbers. Consider the following change of variables: 

u ^ w — hlnu, 

where e R+ \ {0}, > 0; thus u = e^^^, w ER. Denote by the additional 
element — oo and by S the extended real hue R U {0}. The above change of 
variables has a natural extension to the whole 5* by 0^(0) — 0; also, we 
denote D,,{1) = = 1. 

Denote by Sh the set S equipped with the two operations (Bh (gener- 
alized addition) and Qh (generalized multiplication) such that Dh is a ho- 
momorphism of to {S,®h,Qh}- This means that Dh{u\ -\- U2) — 

Dh{ui)®hDh{u2) and D/j(mi-M2) = Dh{ui)QhDh{u2), i.e., wiQhW2 = W1+W2 
and wi (Bh W2 = hhi{e^'^/^ + e^'^^^). It is easy to prove that wi ®h W2 — 
max{wi,W2} as /i — >• 0. 

Denote by Kmax the set 5* = M U {0} equipped with operations © = max 
and = +, where = — 00, 1 = as above. Algebraic structures in M_|_ and 
Sh are isomorphic; therefore Rmax is a result of a deformation of the structure 
in R+. 

We stress the obvious analogy with the quantization procedure, where h 
is the analog of the Planck constant. In these terms, M+ (or M) plays the 
part of a 'quantum object' while Mmax acts as a 'classical' or 'semi-classical' 
object that arises as the result of a dequantization of this quantum object. 

Likewise, denote by Rmin the set R U {0} equipped with operations © = 
min and = -|-, where = +00 and 1 = 0. Clearly, the corresponding 
dequantization procedure is generated by the change of variables u ^ w ~ 
—hlau. 

Consider also the set KU {0, 1}, where = —00, 1 = +cxo, together with 
the operations © = max and = min. Obviously, it can be obtained as a 
result of a 'second dequantization' of M or M+. 

The algebras presented in this section are the most important examples of 
idempotent semirings, the central algebraic structure of Idempotent Analysis. 

Consider a set S equipped with two algebraic operations: addition © and 
multiplication 0. The triple {S, ©, 0} is an idempotent semiring \l it satisfies 
the following conditions (here and below, the symbol -k denotes any of the 
two operations ©, 0): 

• the addition © and the multiplication are associative: x-k {y-k z) = 
{x'ky)i< z for all x,y,z E S; 
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• the addition © is commutative: x (B y = y (B x hi all x, y G S; 



the addition © is idempotent: X © X — X for all X E S; 



• the multiplication is distributive with respect to the addition ©: 

xQ {y® z) = (xQy)® {xQ z) and {x®y) Q z = {xQ z) ®{yQ z) for 
all x,y,z E S. 

A unity of a semiring S is an element 1 G S* such that for all x G S* 

l(Dx = x(Dl = x. 

A zero of a semiring S is an element E S such that 7^ 1 and for all 

X E S 

0®x = x, 0&x = xQO = 0. 

A semiring S is said to be commutative if x & y = y Q x for all x,y E S. 
A commutative semiring is called a semifield if every nonzero element of this 
semiring is invertible. It is clear that Mmax and Rmin are semifields. 

Note that different versions of this axiomatics are used, see, e.g., IQ, 



121 , [T^-[0, IjT^, |2T|], |3y] and some literature indicated in 



these books and papers. Many nontrivial examples of idempotent semirings 
can be found, e.g., in §, §, 0, [g, 0, 0, |2l, 0, 



For example, every vector lattice or ordered group can be treated as an 
idempotent semifield. 

The addition © defines the following standard partial order on S: x ^ y 
if and only ii x(By = y. US contains a zero 0, then ^ x for all x E S. The 
operations © and are consistent with this order in the following sense: if 
X ^ y, then x -k z ^ y -k z and z * x ^ z -k y ior all x,y, z E S. 

The basic object of the traditional calculus is a function defined on some 
set X and taking its values in the field E (or C); its idempotent analog 
is a map X S, where X is some set and 5" = Kmin, Kmax, or another 
idempotent semiring. Let us show that redefinition of basic constructions of 
traditional calculus in terms of Idempotent Mathematics can yield interesting 
and nontrivial results (see, e.g., ||21|, [Q, for details, additional 
examples and generalizations). 

Example 1. Integration and measures. To define an idempotent 
analog of the Riemann integral, consider a Riemann sum for a function if{x), 
X E X = [a,b], and substitute semiring operations © and for operations 
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+ and ■ (usual addition and multiplication) in its expression (for the sake of 
being definite, consider the semiring Mmax): 

AT N 

y](/?(xi)-Ai 1-^ ^(p{xi) Q Ai= max ((p{xi) + Ai), 

where a = < Xi < ■ ■ ■ < xa? = 6, Aj = — Xj-i, i — \^ . . . . As 
maxj Aj ^ 0, the integral sum tends to 



<^(a;) dx — sup ^pix) 

for any function ip: X ^ Rmax that is bounded. In general, for any set X the 
set function 

m^{B) = suY) tp{x) , B C X, 
xeB 

is called an M.„^i^x-i^^(isure on X; since m^pdj^ B^) = snp^m^{Ba)- this mea- 
sure is completely additive. An idempotent integral with respect to this 
measure is defined as 

ip{x)dmu3— / ip{x) Q ^{x) dx — svlY> {'>P{x) + ^{x)) . 

Jx x&X 

Using the standard partial order it is possible to generalize these defini- 
tions for the case of arbitrary idempotent semirings. 

Example 2. Fourier-Legendre transform. Consider the topolog- 
ical group G = M". The usual Fourier-Laplace transform is defined as 

ip{x) ^ (p{i) = I e'^-'^ipix) dx, 

where exp(?;^-.'E) is a character of the group G, i.e., a solution of the following 
functional equation: 

f{x + y)^f{x)f{y). 
The idempotent analog of this equation is 

f{x + y) = f{x)Qfiy) = f{x) + fiy). 

Hence 'idempotent characters' of the group G are linear functions of the form 
X h->- ^ ■ X = ^iXi -|- • • • -|- ^n^n- Thus the Fourier-Laplace transform turns into 

/•e 

(p{x) ^ (p{^) = C-xQ (p{x) dx = sup{^-x + ip{x)). 

JG xeG 
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This is the well-known Legendre (or Fenchel) transform. 

These examples suggest the following formulation of the idempotent cor- 
respondence principle pi]] : 



There exists a heuristic correspondence between interesting, use- 
ful, and important constructions and results over the field of real 
(or complex) numbers and similar constructions and results over 
idempotent semirings in the spirit of N. Bohr's correspondence 
principle in Quantum Mechanics. 

So Idempotent Mathematics can be treated as a 'classical shadow (or 
counterpart)' of the traditional Mathematics over fields. 

In particular, an idempotent version of Interval Analysis can be con- 
structed The idempotent interval arithmetics appear to be remarkably 



simpler than their traditional analog. For example, in the traditional interval 
arithmetic multiplication of intervals is not distributive with respect to in- 
terval addition, while idempotent interval arithmetics conserve distributivity. 
Idempotent interval arithmetics are useful for reliable computing. 



5 Idempotent linearity 

Let be a commutative idempotent semiring. 

The following example of a noncommutative idempotent semiring is very 
important. 

Example 3. Let Matn(S) be a set of all S'-valued matrices, i.e. coef- 
ficients of these matrices are elements of S. Define the sum © of matrices 
A = \\aij\\, B = \\bij\\ G Matn(S) asAQ)B = \\aijQ)bij\\ G Matn(S). The prod- 
uct of two matrices A G Matn(S) and B G Matn(S) is a matrix AB G Matn(S) 
such that AB = \\ aik bkj\\- The set Matn(S) of square matrices is an 
idempotent semiring with respect to these operations. If is the zero of 5*, 
then the matrix O = \\oij\\, where Oij = 0, is the zero of Matn(S); if 1 is the 
unity of S, then the matrix E = \\6ij\\, where dij = 1 if i = j and 6ij = 
otherwise, is the unity of Matn(S). 

Now we discuss an idempotent analog of a linear space. A set V is called a 
semimodule over 5* (or an S'-semimodule) if it is equipped with an idempotent 
commutative associative addition operation (By and a multiplication Qy: S x 
V ^ V satisfying the following conditions: for all X, ^ & S, v, w E V 
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• {X & fi) &v V = X Qv ifJ- &v v); 

• XQv {v ®v w) = (A ©V v) ®v (A Qv w); 

• {X®n)Qvv = {X Qv v) ®v (/^ Qv v). 

An S'-semimodule V is called a semimodule with zero if G and there 
exists a zero element Oy & V such that for all v E V, X E S 

• Ov ®v V = v; 

• XQv Ov = OQv V = Oy. 

Example 4. Finitely generated free semimodule. The sim- 
plest S-semimodule is the direct product 5"" = { (ai, . . . , a„) | aj G S,j = 
1 , . . . , n } . The set of all endomorphisms 5*" —>■ S"^ coincides with the semiring 
Matn(S) of all S"- valued matrices (see example 3). 

The theory of S'-valued matrices is similar to the well-known Perron-Fro- 
benius theory of nonnegative matrices, well advanced and has very many 



applications, see, e.g., §,1, i-|Tl, 0, P, H, [H, 



Example 5. Function spaces. An idempotent function space J^{X; S) 
consists of functional defined on a set X and taking their values in an idem- 
potent semiring S. It is a subset of the set of all maps X ^ S such that if 
fix), g{x) G J^iX; S) and cE S, then (/ © g){x) = f{x) © g{x) E ^(X; S) 
and {cQ f){x) = cQ f{x) E J-'{X; S); in other words, an idempotent function 
space is another example of an S'-semimodule. If the semiring S contains a 
zero element and J-'{X; S) contains the zero constant function o{x) = 0, 
then the function space J^{X; S) has the structure of a semimodule with zero 
o{x) over the semiring S. If the set X is finite we get the previous example. 

Recall that the idempotent addition defines a standard partial order in S. 
An important example of an idempotent functional space is the space B{X; S) 
of all functions X S bounded from above with respect to the partial order 
^ in S. There are many interesting spaces of this type including C{X] S) (a 
space of continuous functions defined on a topological space X), analogs of 
the Sobolev spaces, etc (see, e.g., [|17|, [^, [Q, [^-0] for details). 

According to the correspondence principle, many important concepts, 
ideas and results can be converted from usual Functional Analysis to Idem- 
potent Analysis. For example, an idempotent scalar product in B{X; S) can 
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be defined by the formula 



{(p,^!) = / ip{x) Q ■iIj{x) dx, 

where the integral is defined as the 'sup' operation (see example 1). 

Example 6. Integral operators. It is natural to construct idem- 
potent analogs of integral operators of the form 

K : ip{y) {Klp){x) = / K{x, y) Lp{y) dy, 

JY 

where Lpiij) is an element of a functional space J-'i{Y;S)^ {Kip){x) belongs 
to a space J?-2(X; S) and K{x, y) is a function X xY S. Such operators 
are linear, i.e. they are homomorphisms of the corresponding functional 
semimodules. If S" = Mmaxj then this definition turns into the formula 

{Kip){x) = sup {K{x, y) + (fiy)). 

y& 

Formulas of this type are standard for optimization problems. 



6 Superposition principle 

In Quantum Mechanics the superposition principle means that the Schrodi- 
nger equation (which is basic for the theory) is linear. Similarly in Idempo- 
tent Mathematics the idempotent superposition principle means that some 
important and basic problems and equations (e.g., optimization problems, 
the Bellman equation and its versions and generalizations, the Hamilton- 
Jacobi equation) nonlinear in the usual sense can be treated as linear over 



appropriate idempotent semirings, see |^-||30|, [0. 

The linearity of the Hamilton- Jacobi equation over Rmin (and Mmax) can 
be deduced from the usual linearity (over C) of the corresponding Schrodinger 
equation by means of the dequantization procedure described above (in Sec- 
tion 4). In this case the parameter h of this dequantization coincides with 
i% , where fi is the Planck constant; so in this case h must take imaginary 
values (because /i > 0; see [^] for details). Of course, this is closely related 



to variational principles of mechanics. 

The situation is similar for the differential Bellman equation, see [IT 
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It is well-known that discrete versions of the Bellman equation can be 
treated as linear over appropriate idempotent semirings. The so-called gen- 
eralized stationary (finite dimensional) Bellman equation has the form 

X = AX®B, 

where X, A, B are matrices with elements from an idempotent semiring and 
the corresponding matrix operations are described in example 3 above; the 
matrices A and B are given (specified) and it is necessary to determine X 
from the equation. 

B.A. Carre used the idempotent linear algebra to show that differ- 
ent optimization problems for finite graphs can be formulated in a unified 
manner and reduced to solving these Bellman equations, i.e., systems of lin- 
ear algebraic equations over idempotent semirings. For example. Bellman's 
method of solving shortest path problems corresponds to a version of the Ja- 
cobi method for solving systems of linear equations, whereas Ford's algorithm 
corresponds to a version of the Gauss-Seidel method. 

7 Correspondence principle for computations 

Of course, the idempotent correspondence principle is valid for algorithms as 
well as for their software and hardware implementations |2^-[^. Thus: 

If we have an important and interesting numerical algorithm, then 
there is a good chance that its semiring analogs are important and 
interesting as well. 

In particular, according to the superposition principle, analogs of linear 
algebra algorithms are especially important. Note that numerical algorithms 
for standard infinite-dimensional linear problems over idempotent semirings 
(i.e., for problems related to idempotent integration, integral operators and 
transformations, the Hamilton- Jacobi and generalized Bellman equations) 
deal with the corresponding finite-dimensional (or finite) 'linear approxima- 
tions'. Nonlinear algorithms often can be approximated by linear ones. Thus 
the idempotent linear algebra is a basis for the idempotent numerical analy- 
sis. 

Moreover, it is well-known that linear algebra algorithms are convenient 
for parallel computations; their idempotent analogs admit parallelization as 
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well. Thus we obtain a systematic way of applying parallel computation to 
optimization problems. 

Basic algorithms of linear algebra (such as inner product of two vectors, 
matrix addition and multiplication, etc.) often do not depend on concrete 
semirings, as well as on the nature of domains containing the elements of 
vectors and matrices. Algorithms to construct the closure A* = 1 (B A (B 

© ■ ■ ■ © A" © ■ ■ ■ = v4" of an idempotent matrix A can be derived 

from standard methods for calculating (1 — For the Gauss- Jordan 

elimination method (via LU-decomposition) this trick was used in [^], and 
the corresponding algorithm is universal and can be applied both to the 
Bellman equation and to computing the inverse of a real (or complex) matrix 
(1 — A). Computation of A"^ can be derived from this universal algorithm 
with some obvious cosmetic transformations. 

Thus it seems reasonable to develop universal algorithms that can deal 
equally well with initial data of different domains sharing the same basic 
structure pTl, p2 . 



8 Correspondence principle for hardware de- 
sign 

A systematic application of the correspondence principle to computer calcu- 
lations leads to a unifying approach to software and hardware design. 

The most important and standard numerical algorithms have many hard- 
ware realizations in the form of technical devices or special processors. These 
devices often can be used as prototypes for new hardware units generated by 
substitution of the usual arithmetic operations for its semiring analogs and 
by addition tools for performing neutral elements and 1 (the latter usually 
is not difficult). Of course, the case of numerical semirings consisting of real 
numbers (maybe except neutral elements) is the most simple and natural 
[^-[^. Note that for semifields (including Mmax and Mmin) the operation of 
division is also defined. 

Good and efficient technical ideas and decisions can be transposed from 
prototypes into new hardware units. Thus the correspondence principle gen- 
erated a regular heuristic method for hardware design. Note that to get a 
patent it is necessary to present the so-called 'invention formula', that is 
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to indicate a prototype for the suggested device and the difference between 
these devices. 

Consider (as a typical example) the most popular and important algo- 
rithm of computing the scalar product of two vectors: 

(X, 2/) = + X2?/2 H l-XnVn- (l) 

The universal version of (1) for any semiring A is obvious: 

{X, y) = (Xi © 2/1) © (Xs ?/2) © • ■ ■ © {Xn © Vn)- (2) 

In the case A = Mmax this formula turns into the following one: 

(x, y) = max{xi + yi, Xa + 2/2, ■ " " , a;^ + 2/„}. (3) 

This calculation is standard for many optimization algorithms, so it is 
useful to construct a hardware unit for computing (3). There are many 
different devices (and patents) for computing (1) and every such device can be 
used as a prototype to construct a new device for computing (3) and even (2). 
Many processors for matrix multiplication and for other algorithms of linear 
algebra are based on computing scalar products and on the corresponding 
'elementary' devices respectively, etc. 

There are some methods to make these new devices more universal than 
their prototypes. There is a modest collection of possible operations for 
standard numerical semirings: max, min, and the usual arithmetic opera- 
tions. So, it is easy to construct programmable hardware processors with 
variable basic operations. Using modern technologies it is possible to con- 
struct cheap special-purpose multi-processor chips implementing examined 
algorithms. The so-called systolic processors are especially convenient for 
this purpose. A systolic array is a 'homogeneous' computing medium con- 
sisting of elementary processors, where the general scheme and processor 
connections are simple and regular. Every elementary processor pumps data 
in and out performing elementary operations in a such way that the corre- 
sponding data flow is kept up in the computing medium; there is an analogy 
with the blood circulation and this is a reason for the term 'systolic', see e.g. 



Of course, hardware implementations for important and popular basic 
algorithms can increase the speed of data processing. 
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9 Correspondence principle for software de- 
sign 

Software implementations for universal semiring algorithms are not so effi- 
cient as hardware ones (with respect to the computation speed) but are much 
more flexible. Program modules can deal with abstract (and variable) opera- 
tions and data types. Concrete values for these operations and data types can 
be defined by the corresponding input data. In this case concrete operations 
and data types are generated by means of additional program modules. For 
programs written in this manner it is convenient to use a special techniques of 
the so-called object oriented (and functional) design, see, e.g., [^, [0. For- 
tunately, powerful tools supporting the object-oriented software design have 
recently appeared including compilers for real and convenient programming 
languages (e.g. C++ and Java). 

There is a project to obtain an implementation of the correspondence 
principle approach to scientific calculations in the form of a powerful soft- 
ware system based on a collection of universal algorithms. This approach 
ensures a working time reduction for programmers and users because of the 
software unification. The arbitrary necessary accuracy and safety of numeric 
calculations can be ensured as well. 

The system contains several levels (including programmer and user levels) 
and many modules. 

Roughly speaking, it is divided into three parts. The first part con- 
tains modules that implement domain modules (finite representations of ba- 
sic mathematical objects). The second part implements universal (invariant) 
calculation methods. The third part contains modules implementing model 
dependent algorithms. These modules may be used in user programs written 
in C++ and Java. 

The following modules and algorithms implementations are in progress: 

Domain modules: 

infinite precision integers; 

rational numbers; 

finite precision rational numbers; 

finite precision complex rational numbers; 

fixed- and fioating-slash rational numbers; 

complex rational numbers; 
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arbitrary precision floating-point real numbers; 
arbitrary precision complex numbers; 
p-adic numbers; 
interval numbers; 

ring of polynomials over different rings; 

idempotent semirings i?(max, min), i?(max, +), i?(min, +), interval idem- 
potent semirings 
and others. 

Algorithms: 

linear algebra; 

numerical integration; 

roots of polynomials; 

spline interpolations and approximations; 

rational and polynomial interpolations and approximations; 

special functions calculation; 

differential equations; 

optimization and optimal control; 

idempotent functional analysis 

and others. 

This software system may be especially useful for designers of algorithms, 
software engineers, students and mathematicians. 
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